Geometric phase in dephasing systems 
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Beyond the quantum Markov approximation, we calculate the geometric phase of a two-level 
system driven by a quantized magnetic field subject to phase dephasing. The phase reduces to the 
standard geometric phase in the weak coupling limit and it involves the phase information of the 
environment in general. In contrast with the geometric phase in dissipative systems, the geometric 
phase acquired by the system can be observed on a long time scale. We also show that with the 
system decohering to its pointer states, the geometric phase factor tends to a sum over the phase 
factors pertaining to the pointer states. 
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Quantum mechanics tell us that physical states are 
equivalent up to a global phase, which in general does 
not contain useful information about the described sys- 
tem and thus can be ignored. This is not the case, how- 
ever, for a system transported round a circuit by varying 
the parameters s = (sx,s-z,...) in its Hamiltonian H(s). 
As Berry showed Q , the phase can have a component of 
geometric origin called geometric phase with important 
observable consequences, such as the Aharonov-Bohm ef- 
fect and the spin-^ particle driven by a rotating mag- 
netic field The geometric phases that only depend 
on the path followed by the system during its evolution, 
have been investigated and tested in a variety of settings 
and have been generalized in several directions 3J. The 
geometric phases are attractive both from a theoretical 
perspective, and from the point of view of possible appli- 
cations, among which geometric quantum computation 
is one of the most importance. 



As realistic systems always interact with their envi- 
ronment, the study on the geometric phase in open sys- 
tems become interesting. Garrison and Wright Q were 
the first to touch on this issue by describing open system 
evolution in terms of a non-Hermitian Hamiltonian. This 
is a pure state analysis, so it did not address the prob- 
lem of geometric phases for mixed states. Toward the 
geometric phase for mixed states in open systems, the 
approaches used involve solving the master equation of 
the system 0, l ul fill. Ilil ll^j , employing a quan tum tra- 
jectory analysis |l llll 1 ^ or Krauss operators and the 
perturbative expansions |l7lll8l |. Some interesting results 
were achieved, briefly summarized as follows: nonhermi- 
tian Hamiltonian lead to a modification of Berry's phase 
d, 0|, stochastically evolving ma gnet ic fields produce 
both energy shift and broadening [lfij. phenomenolog- 
ical weakly dissipative Liouvillians alter Berry's phase 
by introducing an imaginary correction or lead to 
damping and mixing of the density matrix elements |l2| . 
However, almost all these studies are performed for dissi- 
pative systems, and thus the representations are applica- 
ble for systems whose energy is not conserved. For open 
systems with conserved energy (dephasing systems), the 



problem beyond the Markov approximation remains un- 
touched to our best knowledge. Because the system- 
environment interaction Hi and the free system Hamil- 
tonian H s commute for dephasing systems, the dynami- 
cal problem and then the geometric phase of the system 
may be solved/calculated precisely. On the other hand, 
the previous study |17| shows that one can not perform 
an arbitrarily long experiment to measure the phase for 
dissipative systems, i.e., it is not allowed to draw phase 
information out of the system on a long time scale, this 
feature of dissipative systems again motivates investiga- 
tion on the problem of geometric phases in dephasing 
systems, where in principle one may get analytical re- 
sults for the phase on any time scale. 

In this Letter, we investigate the behavior of the geo- 
metric phase of a two-level system interacting with a driv- 
ing magnetic field when this system is not only quantized 
but also subjected to decoherence. The environment that 
leads to decoherence may originate from the fluctuation 
in the driving fields, or from the vacuum fluctuations, 
or from the background radiations. We calculate and 
analyze the effect of dephasing of the driven system on 
the geometric phase of the system, our discussions will 
distinguish between two kinds of evolution: (l)The en- 
vironment undergoes an adiabatic evolution, and (2)it 
evolves as it may. 

Let us consider a two-level system driven by a quan- 
tized magnetic field and subjected to decoherence. The 
decoherence process is described by the coupling of the 
two-level system to an environment of harmonic oscilla- 



tors with frequencies 
such a system reads 
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where a, are boson operators for the driving field, er+, 
<7_, and a z are pauli operators for two relevant internal 
atomic levels |e) and \g), and &j, fej are the creation and 
annihilation operators of the environment bosons. The 
system hamiltonian H s characterizes Jaynes-Cummings 
dynamics without neither energy relaxation nor phase 
dephasing, while terms H e and Hi describe a bosonic 
environment and its coupling to the Jaynes-Cummings 
system. The choice of the coupling between the system 
and the environment determines effects of the environ- 
ment. For example, the choice of the system operators 
that do not change the quantum number of the driving 
field when they operate on the dressed state would re- 
sult in relaxations within the dressed states indicated by 
the quantum number n, but not energy relaxations be- 
tween states with different n. The system-environment 
coupling chosen in our model is exactly such a choice, 
and thus, we may rewrite the Hamiltonian Eq. in 
terms of the dressed states as 

h = 0fr n , 

n 

H n = E+(n)\ + (n))(+(n)| 

+ E-{n)\-{n))(-{n)\+Y,^h 



+ vV + I(| + (n))(+(n)| 

- |-M)<-(n)|)5>A J (&}+&,), (2) 



where E±{n) = j" UQ ± hg^/n + 1, and the dressed 
states with indication n are 



to consider the universe (system+environment) to un- 
dergo an adiabatic evolution, in this case the states of 
the environment never evolve during the evolution. Be- 
cause the relaxation of the system due to its coupling 
to the environment is independent of the dressed state 
indication n, the acquired geometric phase of the open 
system are the same as in the case without the environ- 
ment. The second situation is more practical, in which 
we assume the environment initially is in its ground state 
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Hamiltonian Eq 



) C B± . and evolve govern by the 
(JIJ. This initial state means that the 
environment is initialized in the vacuum of modes {bj}, 

or in coherent state | ± Aj ^" +1 )g ± . of modes {B±j}. 
The geometric phase of the universe in this situation may 
be calculated by removing the accumulation of these dy- 
namical phase from the total phase, i.e., 



ff (n) = arg(¥(O)|*Cr)> + i 



( dmt^mt)), (5) 



it is easy to demonstrated that for closed systems <fig(n) 
reduces to the Aharonov-Anandan formula for cyclic evo- 
lutions [2(| and to the Berry phase for adiabatic and 
cyclic evolutions The initial state together with the 
time evolution operator determine the path followed by 
the system, in this sense the geometric phase might de- 
pend on the initial condition. If wc choose 1^(0)) = 
I + W) ® Uj \^)b +>} ( A i = AjVn + l) as the initial 
state, at time t the universe evolves to 



i*(*)> =i+(»)> ®n i— 



(6) 



\±(n)) = -j=(\g,n + l)±\e,n)). (3) 

Here {\n)} stand for the Fock states of the driving field. 
The Hamiltonian H n describes a driven harmonic oscil- 
lator, the driving terms depend on the dressed states, 
but they are independent of the dressed state indication 
n. With these properties, we may simplify the problem 
and restrict our study within the dressed states with the 
same indication n. The eigenstates of the Hamiltonian 
H n take the following form, 



|e-(n)) {n . } 



(»)) ® n 



(4) 



where \nj)B±i denotes the Fock state for new environ- 
ment mode B±j = (bj ± Xjy/n + l/u>j)\ ± (n))(±(n)|, 



i.e., B\ -B 
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>b ± i , depending on the 



dressed states. The corresponding eigenenergies are 
fig^/n + 1 + J2j hujUj and —Tig^Jn + 1 + ■ huijUj, re- 
spectively. We distinguish between two situations to 
study the geometric phase of the system, the first is 



In contrast with classically driving field, in this study 
the driving field is quantized. In order to generate a 
phase change in the state of the field, we borrow the 
idea in Ref.[2l| to introduce the phase shift operator 
U(ip) = exp(—ii(ja'a) and adiabatically apply it to the 
Hamiltonian of the system. Changing ip = Q ■ t slowly 
from to 2n (the corresponding time from to T = ^y) 
the geometric phase generated is calculated by Eq. JSj 
as follows, 
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m=0 



7+(«), 



(7) 



where 7+(n) = (2n + 1)tt is the Berry phase acquired 
when the universe remains unchanged. For a continuous 
spectrum of the environmental modes, the sum over j in 
the above expressions is replaced by an integral involving 
the spectral density with a cutoff frequency ui c 



s{ — ) 2 ,0<lu<uj c , 



(8) 
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this spectrum density is of Ohmic type. Eq. (J3J) and 
Eq.(|HJ together yield 

i _I_ y \ , v W^( n+1 ) T (71 +1) 

= 7+(n)+e c 2 + ^£[coskT) - 1]. 

(9) 

If we choose the eigenstate 1^(0)) = | — (n)) ® Ylj I — 
^-)b_ from another set of eigenstates Eq. J2J as the 
initial condition, the time evolution of the universe can 
be expressed as, 

i 3 

In the same way, we can get the geometric phase pertain- 
ing to this evolution loop, 

<t>g (n) = <f>t( n ) -7+(n) +7_(n), (11) 

namely, the contributions from the system-environment 
coupling are the same for the both pathes. Here 7_ (n) 
is the Berry phase attaining to the dressed state | — (n)) 
in the case without environment, it is easy to prove that 
it takes the same expression as 7+(n). 

The last two terms in equation Eq. @ result from 
the system-environment couplings, they vanish when the 
couplings tend to zero (in the equations, lo c — > 0), thus 
the expressions return to the geometric phases presented 
in Ref. |2l|. The path the system followed changes the 
system-enviroment coupling, and the environment is im- 
possible to return to its initial state due to its huge vari- 
ety of freedom, the geometric phase then depends on the 
time T when we draw out phase information from the sys- 
tem. In the classical limit n — > oo, the contributions from 
the system-environment coupling tend to infinity caused 
by relative strong coupling between the system and the 
driving field, hence it becomes undefined in the sense of 
interferometry. It is interesting to note that for n = 



Here, <& g (n) is the geometric phase for the system under- 
going non- unitary evolution j^. The geometric phase 
Eq. (|13[1 is gauge invariant and can be reduced to the 
well-known results in the unitary evolution, thus it is ex- 
perimentally testable. Now, we exploit the expression 
to calculate the geometric phase for the driven two-level 
system. To this aim, we first write down the reduced 



the phases are not zero, which means that the vacuum 
driving field introduces a correction in geometric phases, 
this expression is relevant when systems are driven by 
fields with few photons. For weak system-environment 
coupling 1^3, Eq. 0( Eq. ltTT|) . in the same way) may 
be expanded in powers of Aj/ujj, up to the second order 
of Aj/uj, Eq.0 follows, 

t g (n)= 1+ (n)+J2^T(^f. (12) 

— LOj 

3 J 

The explanation of this result is very simple, as the 
system-environment coupling is very weak, the most con- 
tribution to the geometric phase come from the system- 
driving field coupling. The same dependence of 4> g (n) on 
Aj (i.e., Ap no contribution proportional to Aj) can be 
found in Ref. |l7| . 

Up to now, we have calculated the geometric phase for 
the universe, explicit expressions for the geometric phase 
were obtained, these expressions are of relevance for the 
case in which systems are coupled to an environment that 
describes parameter fluctuations. For instance, imperfect 
dipole transitions between states \g) and |e) due to fluc- 
tuation of the driving laser intensity may be modelled 
by the coupling in Eq. in this situation the environ- 
ment and the driving field are the same. This is not the 
case, however, when the environment is an independent 
system (say, black body radiations), we have to trace 
out the environment in order to calculate the dynami- 
cal information for the system, thus the evolution of the 
system is no longer unitary. For non-unitary evolution, 
the geometric phase can be calculated as follows. First, 
solve the eigenvalue problem for the reduced density ma- 
trix p(t) and obtain its eigenvalues £fc(i) as well as the 
corresponding eigenvectors \ipk(t)); secondly, substitute 
e k (t) and |^fc(*)) into 



(13) 

I 

density matrix p(t), 

*>-U£ w c tT> ( ")- <"> 

where the initial state of (c+| + (n)) + c_| — (n))) <E> 
Yij i s assumed for the universe, and F(t) = 

exp[-j-(E + (n) - E_{n))t] ■ exp[- £j i&(t)] with Vj (t) = 
4|^j-| 2 (l — cos ujjt). Simple algebra gives the eigenvalues 
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and the corresponding eigenvectors of p(t), 

e±(t) = \±y{\c + \*-\c^y+±\c%^F(t)\\ 
\i>±{t)) = X ± {t)\ + {n))+Y ± {t)\ - (n)>, (15) 

where X ± (t) = c* + c-F(t)/y/ \d* + c-F(t)\* + (e ± (t) - \c+\*), 
and Y ± (t) = ^1 - \X ± {t)\ 2 . Eq. JTHJ) and Eq.lO to- 
gether yield the geometric phase for the system. The 
expression for the geometric phase is tedious, so instead 
of writing down the expression, we present here some 
remarkable comments. The expressions for the geometric 
phase are analytically exact, so we might predict the 
behavior of the geometric phase on a long time scale. For 
any j, r)k(t) > 0, so with t — > oo, F(t) — * for a random 
spectrum density of the environment, this indicates that 
$ 9 (T -> oo) -> arg[|c+| 2 e-^+M + jc-^e-^-W] for a 
relative long T. This result is quite interesting: With 
the off-diagonal elements of the reduced density matrix 
tending to zero (decoherence), the driven two-level 
system decohers to its pointer states | + (n)) or | — (n)), 
while the geometric phase factor of the driven system 
reduces to a weighted sum over the phase factors 
pertaining to the pointer states, this provides us a new 
way to observe decoherence effects. 

Summarizing, the geometric phases for a dephasing 
system (open system) have presented and discussed. 
The open system has been demonstrated by a driven 
two-level system coupling to an environment of har- 
monic oscillators. The results show that there are no 
correction in the geometric phase of the universe due to 
the system-environment coupling when the environment 
undergoes an adiabatic evolution, whereas the correction 
in the phase is path dependent when the constraint on 
the evolution is released. The mixed state geometric 
phase for the dephasing system is also presented and 
discussed. The geometric phase factor would tend to a 
sum over these phase factor pertaining to the pointer 
states with the open system decohering to its pointer 
states, it is a reflection of the decoherence in geometric 
phases. 
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